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We calculate the large-g series expansion for the ground state degeneracy 
exponent of the ground state entropy) per site of the g-state Potts antifer- 



romagnet on the (3 • 12 2 ) lattice, to order 0(y 19 ), where y = l/(q — 1). We 

s. 

note a remarkable agreement, to 0(y 18 ), between this series and a rigorous 



lower bound derived recently. 
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I. INTRODUCTION 

Nonzero ground state disorder and associated entropy, So 7^ 0, is an important subject in 
statistical mechanics. One physical example is provided by ice, for which the residual molar 
entropy is S = 0.82 ±0.05 cal/(K-mole), i.e., S /R = 0.41 ±0.03, where R = N Avog k B [|,|. 
A particularly simple model exhibiting ground state entropy without the complication of 
frustration is the g-state Potts antiferromagnet (AF) H[| on a lattice A, for q > x(A), where 
x(A) denotes the minimum number of colors necessary to color the vertices of the lattice such 
that no two adjacent vertices have the same color. This model has a deep connection with 
graph theory in mathematics, since the zero-temperature partition function of the above- 
mentioned g-state Potts antiferromagnet on a lattice A satisfies Z(A, q,T — 0)paf = P{A, q), 
where P(G, q) is the chromatic polynomial |j| expressing the number of ways of coloring the 
vertices of a graph G with q colors such that no two adjacent vertices (connected by a bond 
of the graph) have the same color. Hence, the ground state entropy per site is given by 
So/ks = In W(A, q), where W(A, q), the ground state degeneracy per site, is 

W(A,q)=\imP(A n ,q) 1 / n (1.1) 

Here, A n denotes an n- vertex lattice of type A with appropriate (e.g., free) boundary condi- 
tions. Since nontrivial exact solutions for this function are known in only a very few cases 
(square lattice for q = 3 ||, triangular lattice f7j, and kagome lattice for q = 3 HU), it 
is important to exploit and extend general approximate methods that can be applied to all 
cases. Such methods include rigorous upper and lower bounds, large-g series expansions, 
and Monte Carlo measurements. Recently, with R. Shrock, the present author studied the 
ground state entropy in antiferromagnetic Potts models on various lattices and obtained 
further results with these three methods ||- [iT|. We derived a general lower bound on 



W(A, q) |TJ] which applies to all Archimedean lattices and coincides to many orders with 
large-g series expansion of this function. Previous large-g series expansions include works 
by Baker |12[|, Nagle [O,0|, Kim and Enting |15||, Bakaev et al [ffB] and work reported in 



I jnjJnj . Related work on series expansions for the ground state degeneracy of ice was done 



by Nagle [17 



Large-g series expansions of the respective W(A, q) functions on various Archimedean 
lattices were computed in Ref. [O]. In particular, W(A,q) for A = (3 ■ 12 2 ) was computed 



to 0(y 13 ) [TT[|. In the present paper we extend this series to higher order, namely to 0(y 19 ). 
Our main motivation is to check the accuracy of the lower bound on W((3 • 12 2 ),g) given 
in |ll[| . It is interesting that this lower bound coincides with the first 19 terms, i.e. to 
0(y 18 ), in the large-g series. We choose the lattice A = (3 • 12 2 ) as an illustrative example of 
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a heteropolygonal Archimedean lattice. The reader is referred to Refs. ||- [|TTJ for further 
background and references. 



II. LARGE-g SERIES EXPANSION 

Before proceeding, we recall that an Archimedean lattice is defined as a uniform tiling 
of the plane by regular polygons in which all vertices are equivalent [pl|. Such a lattice 



is specified by the ordered sequence of polygons that one traverses in making a complete 
circuit around a vertex in a given (say counterclockwise) direction. This is incorporated in 
the mathematical notation for an Archimedean lattice A: 

A = (IK*) (2.1) 

i 

where in the above circuit, the notation p"* indicates that the regular polygon pi occurs 
contiguously ai times; it can also occur noncontiguously. We shall denote a^ s as the sum of 
the Oj's over all of the occurrences of the given pi in the product. Because the starting point 
is irrelevant, the symbol is invariant under cyclic permutations. The number of polygons of 
type pi per site is given by 

u Pi = ^ (2.2) 

Pi 

The coordination number for an Archimedean lattice is A = J2% a i,s- I n particular, for the 
(3 • 12 2 ) lattice considered in this paper, the number of triangles per site is p% = 1/3, the 
number of 12-gons per site is pu = 1/6, and the coordination number is A = 3. A section 
of this lattice is shown in Fig. ([!]) . 

A general upper bound on a chromatic polynomial for an n- vertex graph G is P(G, q) < 
q n . This yields the corresponding upper bound W({G},q) < q. Hence, as in our previous 
work ||- |lTJ, it is natural to define a reduced function that has a finite limit as q — ► oo, 



W r ({G},q) = q- l W({G},q) (2.3) 

When calculating large-g Taylor series expansions for W functions on regular lattices , it is 
most convenient to carry this out for the related function 

for which the large-g series can be written in the form 



with 



W(A, y ) = i+J2 WK, m y r 



m=l 



(2.5) 



q-1 



(2.6) 




FIG. 1. Section of the (3 • 12 2 ) Archimedean lattice 



Our calculations of large-g series expansion use the method of Ref . |TJ| . The chromatic 
polynomial is written as the sum 



P(G,q) 



{q ~ 1)E E (-i) e q(e ~ v) 



JE-n) 



G a <G 



(9-1) 



■m(G a ,q), 



(2.7) 



where E is the total number of edges of the n-vertex graph G, m(G a , q) are weights [[0| of 
weak subgraphs G a of G, and e and v are the numbers of edges and vertices, respectively, of 
G a . The summation is over all weak subgraphs G a . The weight function m(G a , q) vanishes 
if G a has any vertices of degree one or if G a has a bridge. Another property is that m(G a , q) 
does not change under the insertion or deletion of vertices of degree two in G a . Thus in the 
summation of weak subgraphs, one effectively has to consider (connected and disconnected) 
subgraphs with no vertices of degree one and without bridges. One also has to consider 
graphs with articulation points. Weight functions m(G a , q) are independent of G and satisfy 



a simple recursion relation formula. Ref. |L4| gives the weights of all star graphs with 
cyclomatic number less or equal to four. p0| 

The series expansion for VF((3 • 12 2 ),g) to 0(y 19 ) involves star graphs with cyclomatic 
number up to 7. The graphs with cyclomatic number equal to 5, 6 and 7 which enter in the 
series expansion to this order are shown in Figs.(@a), (|]d) and (|]e), respectively. To derive 



the weights of these graphs, we use Theorems III and VII of Ref. ||14jl . Theorem III states 
that if a graph G consists of two pieces G\ and G2 which have just one vertex in common, 
its weight is given by 

m(G, q) = -m(Gi, q)m(G 2 , q). (2.8) 

Q 

Theorem VII states that 

m(G,q) = --m(G',q) + m(G",q), (2.9) 

Q 

where G' is derived from G by omitting the edge between two vertices, say i and j, and G" 
is the graph with vertices i and j identified. As an example, consider the graphs G5, G' 5 and 
G'l depicted in Figs.(|a), (|2|b) and (|c). G'l has an articulation point and, using Eq. (|2.8|) , 
we can write its weight as m(G'l, q) = q~ l m(P, q)m(G' 5 , q), where P here stands for polygon. 
The weights for P and G' 5 are m(P,q) = (q — 1) and m(G' 5 ,q) — (g — l)(g — 2) 3 /g 3 [|14 |. 
Hence, Eq.( [2.9|) yields 



m(G 5 , q) = --m(G' 5 , q) + m{<%, q) = \{q - l)(g - 2) 4 . (2.10) 

Note that vertices of degree two have been omitted in Fig.(Q). 

The weights of the graphs with higher cyclomatic numbers, shown in Figs.(^d) and (0e), 
can be similarly determined to be 

m(G 6 ,q) = ^(q-l)(q-2) 5 , (2.11) 

q° 



and 



respectively. 



m(G 7 ,q) = -(q-l)(q-2) 



(2.12) 



i J 




(a) G, 




(b) G' 




(c) G 





(d) G, 



(e) G 



FIG. 2. Illustration of star graphs with cyclomatic number c < 7. (a) G5 is a graph with c = 5; 
(b) G' 5 is derived from G5 by omitting one edge between vertices i and j, and (c) G'l is the graph 
with vertices i and j identified, (d) G% is a graph with c = 6 and (e) Gj has c = 7. Vertices of 
degree two are not shown. 






(a) 



(b) 



(c) 






(d) 



(e) 



(f) 



FIG. 3. Illustration of some graphs which enter in the large-g series for W((3 • 12 2 ),q). Graphs 
(a), (b) and (c) enter in series to 0(y 12 ), while graphs (d), (e) and (f) contribute to 0(y 13 ) and 
higher. 

The subgraphs which contribute to the series to 0(y 12 ) are (i) graphs formed by t discon- 
nected triangles, where t = 1,2, ...,6, (ii) polygons with 12 vertices (12-gons), as shown in 
Fig.(^a), (iii) graphs formed by a triangle connected to a 12-gon, as shown in Figs.(|^b) and 
(0c). Let us refer to the edge in the overlap between a triangle and a 12-gon as an internal 
edge. Internal edges can be part of the subgraph (as in Fig.(^c)) or not (as in Fig. (0b)). 

Further subgraphs which contribute to 0(y 13 ) are (i) graphs formed by one 12-gon and 
one disconnected triangle and (ii) graphs formed by two triangles connected to a 12-gon. 
In the latter case, one has to consider graphs with no internal edges (as in Fig.(|3]d)), one 
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internal edge (as in Fig.(]3|e)) and two internal edges (as in Fig.(^f)). Moreover, one has to 
consider all the distinguishable permutations in the positions of the two triangles. In the 
remaining of this paper, when we refer to t triangles connected to 12-gons, we are including 
all possible distinguishable permutations of the triangles and all cases where i internal edges 
belong to the subgraph, with i — 0, 1, ..., t. 

To 0(y 1A ) other subgraphs which enter in the series are (i) graphs formed by 7 discon- 
nected triangles, (ii) graphs formed by 3 triangles connected to a 12-gon, (iii) graphs formed 
by one triangle connected to a 12-gon and one disconnected triangle. 

Series to 0(y 15 ) includes (i) graphs formed by 4 triangles connected to a 12-gon, (ii) 
graphs formed by 2 triangles connected to a 12-gon and one disconnected triangle, (iii) two 
disconnected triangles and one 12-gon. 

To 0(y m ), extra subgraphs are (i) graphs comprised of 8 disconnected triangles, (ii) 
graphs formed by 5 triangles connected to a 12-gon, (iii) graphs composed of a triangle 
connected to a 12-gon and two disconnected triangles, (iv) graphs formed by 3 triangles 
connected to a 12-gon and one disconnected triangle. 

To 0(y 17 ), other subgraphs are (i) graphs formed by 6 triangles connected to a 12-gon, 
(ii) 3 disconnected triangles and one 12-gon, (iii) 2 triangles connected to a 12-gon and two 
disconnected triangles, (iv) 4 triangles connected to a 12-gon and one disconnected triangle. 

To 0(y 18 ), further subgraphs are (i) graphs formed by 9 disconnected triangles, (ii) 
graphs formed by one triangle connected to a 12-gon and 3 disconnected triangles, (iii) 
graphs comprised of 3 triangles connected to a 12-gon and two disconnected triangles, (iv) 
graphs formed by 5 triangles connected to a 12-gon and one disconnected triangle. 

To 0(y 19 ), new subgraphs are (i) graphs formed by 4 disconnected triangles and one 12- 
gon, (ii) graphs formed by 2 triangles connected to a 12-gon and 3 disconnected triangles, (iii) 
graphs formed by 4 triangles connected to a 12-gon and 2 disconnected triangles, (iv) graphs 
formed by 6 triangles connected to a 12-gon and one disconnected triangle, (v) 20-gons. 

To this order, we obtain 

W((3 ■ 12"),,,) = 1 - \y> - \y> - Lf - f / - §,» + i„" - %y» 

- r/ 3 - fy" - h^ - f^ 16 - l e y" ~ ^ 18 + 1^ 19 + °^ < 2 - 13 > 



The lower bound of Ref . |T]J , namely 

W((3 • 12 2 ),2/), = (1 - y 2 ) 1 /^! + y 11 ) 1 / 6 , (2.14) 



coincides with the first 19 terms of the series given in Eq. (|2.13|) , i.e. to 0(y ls ). This is 



remarkable and shows that this lower bound is indeed a very accurate approximation to the 



exact solution for the W function. The lower bound first differs from the large-g expansion 
for the exact W function at order y 19 : the Taylor series expansion of this lower bound gives 
— (5/729)?/ 19 whereas the large-g series expansion of W yields (719/1458)?/ 19 . 

It is interesting to note that the lowest order in y in which the bound ( p.!4j ) differs from 
the series is an order in which subgraphs involving two adjacent 12-gons, i.e. 20-gons, first 
contribute in the series expansion. If one were to calculate the series expansion without 
considering the contribution of 20-gons to 0(y 19 ), one would get a result which coincides 
with the coefficient of the 0(y 19 ) term in the Taylor series of the bound ( 2.14Q . 



III. CONCLUSIONS 

We report on large-g series expansion for the ground state degeneracy of the Potts anti- 
ferromagnet on the (3- 12 2 ) lattice, to 0(y 19 ). It is remarkable that the lower bound derived 
previously coincides with the first 19 terms of the series, i.e. to 0(y 18 ). It is worthwhile to 
perform similar series expansions to high orders for other types of lattices. 

I would like to thank Professor Robert Shrock for helpful comments on the manuscript. 
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